Resonant X-ray scattering investigation of the multipolar ordering in Ca 3 Co206 
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We have performed a resonant x-ray scattering (RXS) study near the Co K edge on a single 
crystal of Ca3Co2C>6. In the magnetically ordered phase a new class of weak reflections appears 
at the magnetic propagation vector f ~ (31515)- These new reflections allow direct access to the 
dipolar-quadrupolar E1E2 scattering channel. The theoretical possibility of observing isolated E1E2 
electromagnetic multipoles has attracted a lot of interest in the recent years. Unfortunately in many 
system of interest, parity even and parity odd tensor contributions occur at the same positions in 
reciprocal space. We demonstrate that in Ca3Co2 0e it is possible to completely separate the parity 
even from the parity odd terms. The possibility of observing such terms even in globally centrosym- 
metric systems using RXS has been investigated theoretically; Ca3Co2 0e allows a symmetry based 
separation of this contribution. 

PACS numbers: 75.25.+Z, 75.50.Ee, 78.70.Ck 
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Ca3Co-206 has until very recently been described in 
terms of ferromagnetic chains coupled antiferromagneti- 
cally on a triangular lattice. The discovery of a small in- 
commensurability [l[ in the magnetic propagation vector 
along the c axis has been crucial in leading to a reeval- 
uation of the role of the local exchange integrals in the 
description of the system. The system, shown in Fig. 1, 
consist of chains made up of alternating distorted octa- 
hedra and trigonal C0O6 prisms sharing faces, running 
along the hexagonal c axis and arranged in a triangu- 
lar pattern within the ab plane Q. The different local 
environments leave the Co 34 " ions on the octahedral site 
(Col) in a low-spin (S = 0) state, and those on trigonal 
prism (Coll) sites in the high-spin (S = 2) state [H, 0). 
The local anisotropy of the trigonal prism forces the mag- 
netic moments to point along the c axis as confirmed by a 
number of experimental results. Ca3Co20g is usually de- 
scribed in the hexagonal setting of the R3c space group, 
as this representation allows one to immediately identify 
the triangular arrangement of the CoOg chains within 
the ab planes. However, as this setting is non primitive, 
it makes both the description of the magnetic structure 
and the analysis of the symmetry of the tensors more 
difficult. For these reasons, in this letter the rhombohc- 
dral axes are used throughout. In the rhombohcdral set- 
ting, the unit-cell dimensions are a=b=c=6.274 A and 
a = [3 = 7 = 92.53°. Compared with other transition 
metal ions, the Co 3+ ion carries a large orbital moment 
that in Ca 3 Co 2 06 reaches 1.7 Hb with L/S ~ 0.94 
This large orbital contribution is crucial to the observa- 
tion of small effects such as the E\Ei signal described in 
the following. 

In this letter we report the results of a resonant x- 
ray investigation on Ca3Co206 performed near the Co 
K-edge. In the magnetically ordered state, in addition 
to the principal magnetic reflections, we observe a sec- 




FIG. 1: (Color on line) Schematic of the unit cell of 
Ca3Co2 0e. The trigonal prisms are shown in light blue (dark 
grey) and the octahedra in yellow (light grey). The magnetic 
moments are also drawn. 



ond class of reflections that also appear at the magnetic 
propagation vector f ~ (3, g, g) but at positions where 
the magnetic structure factor is zero. These reflections 
are characterized by a completely different photon energy 
spectra and azimuthal dependence. We first demonstrate 
that unlike the principal magnetic reflections, they can- 
not be described within the dipolar-dipolar (EiE\) and 
quadrupolar-quadrupolar {E2E2) approximation. We 
then provide a fit based on a magnetic contribution com- 
ing from the dipolar-quadrupolar (E1E2) interference. 
The occurrence of higher rank scattering in the space 
group i?3 c was treated by a number of authors @, 0, S @] • 
The case of V2O3 [loj [. where a new family of reflec- 
tions was discovered, looks similar to our situation. A 
very important difference between Ca3Co206 and V2O3 
is the presence of a first order phase transition in V2O3. 
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In fact, there is no symmetry left at the V site in the 
low temperature phase of V2O3. For Ca 3 Co206, no de- 
tectable structural changes are observed as the system is 
cooled, even though the magnetic symmetry of the sys- 
tem is lower than the structural symmetry. Therefore, 
we preserve all the symmetry of the space group when 
we consider the charge distribution contribution to the 
scattering. 

Single crystals of Ca 3 Co20g were grown by a flux 
method. Two crystals with different orientations were 
used for the RXS experiments. The first crystal was 
5x2x1 mm 3 with the largest natural face per- 
pendicular to the (1 1) reflection. The second crystal 
was polished to have a diffraction face perpendicular to 
the (1,1,1) direction. The RXS experiments were per- 
formed at the magnetic scattering beamline ID20 [11| at 
the ESRF (Grenoble). The beamline optics were opti- 
mized at 7.7 keV, close to the Co K-edge. The sample 
was mounted in a displex cryostat. The diffractomctcr 
was operated in the vertical plane scattering mode with 
an azimuth set-up to allow for a sample rotation about 
the scattering vector. Hence, the natural polarization 
of the incident beam was perpendicular to the scatter- 
ing plane (a) . The integrated intensity of the reflections 
was measured using a photon counting avalanche photo- 
diode detector. The polarization of the reflected beam 
was linearly analyzed by rotating the scattering plane of 
a highly oriented (00L) pyrolitic graphite plate. 

In the resonant regime, the elastic scattering ampli- 
tude, A(Q,u>), can be written in terms of the tensorial 
atomic scattering factor fj(u>) with the index j running 
over all the atoms contributing to the scattering in the 
unit cell: 
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where TiQ is the momentum transfer in the scattering pro- 
cess and hut is the energy of the incoming and outgoing 
photons. The two magnetic Co ions in the unit cell lie on 
the trigonal axis at the special positions Co\ = 
and Co 2 = (1,1,1), related by the inversion centre / at 
(0, 0, 0). We note that due to the small incommensurabil- 
ity of the magnetic propagation vector, the value of the 
moment at site two is not identical to the moment at site 
one, therefore the inversion symmetry is locally broken. 
However, as the difference between the moment on the 
two ions is small, we neglect this difference and we keep 
this symmetry operation. As the magnetic structure of 
the system is incommensurate, it is still possible to find 
one inversion point in the magnetized crystal. Therefore 
globally the system cannot be ferroelectric. 

The non-magnetic Co ions that may contribute to a 
non-magnetic resonant process lie on the inversion center 
(Co™" = (0, 0, 0), Co™" = {\,\, 5)) and they are related 



by an axis 2. In the following, we neglect these ions as 
an analysis of the tensors contributing to the scattering 
shows that the only tensors that may contribute to the 
scattering at the reflection positions (F 3 + F^ 3 ) produce 
an azimuthal dependence that does not match with the 
experimental observations. 

The scattering amplitude due to the magnetic Co ions 
is given by: 

A(Q,lo) = e^-Ciw) (/ x + J/ ie **<*'*.*)) . (2) 

The tensorial atomic scattering factor fj (u>) is a complex 
quantity whose explicit form is given in many papers. 
Following Di Matteo et al. Q it is useful to write 

f j (cj) = J2(-) q TVF* q (j;u), (3) 
p. 9. 

where the irreducible tensors describe the scattering 
process and are a function of the incident and outgoing 
polarizations and wave vectors, while the irreducible ten- 
sors of the same rank F^_ q (j;co) describe the multipolar 
properties of the system. It is useful to evaluate the ten- 
sors and F p _ q (j;co) in two different reference frames 
and to rotate F^_ q (j; to) to bring the tensors into the same 
reference frame. are calculated in the reference frame 
given in Ref. @, whereas the crystal reference frame has 
been chosen with the z axes parallel to the (111) direc- 
tion, the axis y parallel to the (101) direction and the 
axis x parallel to y x z. 

The F^_ tensors allowed by symmetry need to be in- 
variant under the point group of the system. The non- 
magnetic point group symmetry is 32. We note that 
maintaining this point group would lead to the absence 
of a Ei Ei magnetic contribution as expected since 32 is 
not an admissible magnetic point group. 

The admissible magnetic point group is 32'. The axis 
3 is maintained and since the proper rotations act on 
spins in the same way as on polar vectors, we need to 
reverse the action of the axis 2 to maintain the spin in- 
variance under the point group. The tensors that survive 
to the axis 3 are as follows: (i) in the E\Ei channel: 
F°,F^,Fg; (ii) in the E X E 2 channel: Fq 1 , # 2 , F 3 , F$ ± 
F* 3 , and F^,FgiF§iF$ ± F 3 3 ; (iii) in the E 2 E 2 chan- 
nel: F °, Fi,Fl Fl F 3 ± F 3 3 , F 4 , F 3 4 ± F± 3 . The tilde 
indicates the time reversal odd (i.e. magnetic) tensors. 

A further reduction in the number of the tensors con- 
tributing to the scattering is obtained by applying the 
axis 2 on the non-magnetic tensors and 2' if we consider 
the magnetic quantities. We are left with: (i) E\E\\ 
Fl FlFl, (ii) EiE 2 : F \ F§, F 3 -F 3 3 and F* 
(iii) E 2 E 2 : F°, Fl Fl + F 3 3 , F 4 , F 3 ~ F 4 3 . It is easy to 
see Fq is retained and that it corresponds to the magnetic 
moment along the trigonal axis. 
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FIG. 2: (Color on line) The intensity versus photon-energy 
dependence around the Co K-absorption edge of an E±Ei 
(top) and an E\Ei (bottom) magnetic reflection. The data 
were collected in the on channel. 



There arc still a considerable number of contributions. 
An effective way to simplify the problem is to calcu- 
late the structure factor in Eq. [5] for the reflections 
Q = (H,K,L) ± (t, t, t), with r = i for the magnetic 
reflections and t = for the structural reflections. In 
the magnetic case we start with H + K + L = even, 
then H + K + L ± 3r = odd and A(hkl ± f, w) = 
±i(Fq 1 — (— l) p i 7 'g 1 ), where P is the parity of the ten- 
sor. Only the parity odd tensor, E\E 2 processes, will 
contribute to the scattering process. Considering the for- 
bidden reflections H+K+L = odd, then H+K+L±3t = 
even and A(hkl±T,u) = ±(F (J p 1 + (-l) p i^ p 1 ), we see that 
only the E^Ey and the E 2 E 2 processes will contribute to 
the scattering process at these positions. 

In Fig. [5] we report the photon-energy dependence 
around the Co K-absorption edge of the (§,§,§) and 

the (|, |, I) reflections. These two reflections are the 
(3,3, 2) — (|, |, |) and the (3,0, -2)-(|, ±, §) reflections 
respectively, so the first one provides us with access to 
the E\Ei terms and the second main magnetic reflection 
allows us to probe the E\E\ and E2E2 contributions. 

The main magnetic reflections (e.g. (|, ^, |)) have a 
complex energy spectra that has been discussed in ref. [l[ . 
The new (|, |, |) reflection exhibits a very simple energy 
spectra that can be described using a single oscillator 
localized in the pre-edge region, which points to the in- 
terference term (E1E2) being well localized in energy in 
the pre-edge region. This is in good agreement with the 
fact that the p states will be not accessible through these 
higher rank processes. The reduced intensity (~ j^) of 

the (§,§,§) reflection compared to the (|, |, |) reflec- 
tion confirms that the former originates from a higher 
rank process. 

Fig. [3] shows the temperature dependence of the main 
magnetic reflection (§,-|,|) and of the EiE 2 reflection 

(if' 5> §)■ ^ can ^ e seen *^at ^ ne re fl ec ti° ns not only 
appear at the same magnetic propagation vector as the 
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FIG. 3: (Color on line) Comparison of the temperature de- 
pendence of two reflections due to different scattering pro- 
cesses. The lines are a guide to the eye. While the general 
behavior is similar, the intensity of the E1E2 reflection (cir- 
cles) increases more slowly immediately below the ordering 
temperature than the signal due to E1E1 (squares). 

main magnetic reflections, but that it also exhibits a very 
similar temperature dependence. These two facts suggest 
that the two families of reflections share the same physi- 
cal origin. It is worth noting that due to the sensitivity 
of anomalous diffraction to small local structural changes 
we cannot completely exclude the possibility that an 
atomic displacement may follow the magnetic ordering. 
However, no displacements have been reported in this 
compound. We therefore assume that only the time odd 
tensors can contribute to the two family of reflections. 
The azimuthal dependence of the reflections reported in 
Fig. gjconfirms that the (§, §, §) and the |, |) reflec- 
tions are due to higher rank tensors. In fact the (|, |, |) 
reflection due to the EiEx process has an almost perfect 
2-fold symmetry, whereas the (§,§,§) signal, which is 
« 10 degrees from the 3-fold (111) axis has 6 peaks with 
different intensities and with zeros values that are not ex- 
actly 60 degrees apart. The term F3 — E^_ 3 has a 3-fold 
periodicity with respect to a rotation about the quantiza- 
tion axis, whereas Fq and F§ are constant therefore they 
cannot reproduce the symmetry of the signal. Hence we 
will neglect these contributions in the data analysis. 

To calculate the behavior of the reflections in the ref- 
erence frame of we apply a sequence of coordinate 
rotations Q to the tensorial structural factor. This op- 
eration, together with the actual calculations of the in 
terms of the initial and final polarization and wave vec- 
tors, produces the following expressions for the (§,§,§) 
signal 

(jo- = a 5 T| [a 6 + a 7 C2^]c0S^ (4) 

(77T = aiT^c^cg [a 2 C6/s^ + (a 3 + a 4 C2^)se] (5) 

where are complex numerical coefficients defined by 
the direction in the space of the reflection, c x (s x ) are 
shorthand for cosa;(sina:), 9 is the Bragg angle and -0 is 
the azimuthal angle. The variable T^(cu) embeds the res- 
onant denominator and the radial part of the integrals 
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FIG. 4: (Color on line) Theoretical and experimental az- 
imuthal dependence of three magnetic reflections at the Co 
K-absorption edge. In the first panel an E±Ei reflection is 
reported, whereas the central and bottom panels show reflec- 
tions resulting from the E\E^ interference together with the 
model predictions. The amplitude is the only free parameter 
in the fit. 



between the ground state and the intermediate state in 
the resonant process. A similar expression is obtained 
in the case of the (^p §> I) reflection. As the system 
is not magneto-electric, in an absorption experiment it 
would not be possible to observe a tensorial contribu- 
tion coming from the E1E2 interference. However, in a 
resonant diffraction experiment, we need to look at the 
local symmetry to establish which tensors can contribute 
to the scattering. The possibility of observing an E1E2 
interference term in a globally ccntrosymmctric system 
has been demonstrated in the case of V2O3 in the time 
even case, where the absence of an inversion symmetry at 
the V position is responsible for the appearance of such 
a term. In CasC^Oe at the Coll site, both the mag- 
netic 32' and the non-magnetic point groups 32 allow a 
magneto-electric tensor to occur. A useful classification 
of the electromagnetic terms occurring in the RXS scat- 
tering based exclusively on the linear magneto-electric 
effect has been given in Refs. 



12L [13. According to this 



scheme, the dominant term contributing to the scatter- 
ing is a polar toroidal octupole (anapole). A contribution 
from a polar toroidal dipolc (Tq) is also allowed. 

The possibility of observing a polar toroidal octupole 
(anapole) has been widely discussed theoretically, but to 
date very few experimental observations of this quantity 
are available. V2O3 is certainly the most well known case 
[nj but in this instance the reflection structure factor 
allows an E2E2 term (the magnetic octupole) to occur 
together with the E1E2 terms 12(. More recently, a pos- 
sible E1E2 signal has been reported in the ferroelectric 



phase of TbMn0 3 [14f. Here, both the E X E 2 and the 
Ei Ei contributions are symmetry allowed, but occur at 
different energies. So the tail of the £4 £4, that is much 
stronger in intensity, can give a significant contribution 
at the lower energies where the E1E2 resonates. 

In Ca 3 Co206 the structure factor allows one to sep- 
arate completely the parity odd from the parity even 
events. The importance of the E1E2 interference in the 
analysis of RXS II 311 and of the anapole in the electro- 
magnetic theory [l5| as well as in nuclear and particle 
physics has been underlined in a number of experimental 
and theoretical papers. Ab-initio simulations have re- 
cently revealed the existence of toroidal polarizations in 
nanodiscs and nanorods[3| and their use was suggested 
for data storage. Usually the existence of such effects 
is masked in solid state physics by other effects. In this 
respect RXS offers an almost unique possibility of access- 
ing such effects. Other sources of scattering (e.g. non- 
magnetic E1E2 interference) cannot be ruled out by the 
geometrical dependence of the RXS signal alone. Nev- 
ertheless, the absence of any structural change and the 
thermal behavior of the signal strongly suggests that a 
form of anapolar ordering develops in the low tempera- 
ture phase of Ca 3 Co 2 06. In conclusion, we demonstrate 
once again that RXS is a powerful technique for access- 
ing subtle features in the charge distribution in complex 
materials and that there is an urgent need to develop 
further the theoretical interpretation of the data. 
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